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Abstract 

We determine running coupling corrections to the kernel of the non-linear evolution 
equation for the cross section of single diffractive dissociation in high energy DIS. The 
running coupling kernel for diffractive evolution is found to be exactly the same as the 
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kernel of the rcBK evolution equation. 
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&■ 1 Introduction 

In recent years running coupling corrections have been calculated for a variety of small- a; 
observables. The running coupling corrections for the linear Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) evolution equation [T}[2] were found in [3] (see also [HE]), while the corrections for 
the non-linear Balitsky-Kovchegov (BK) [6rfT0] and Jalilian-Marian-Iancu-McLerran-Weigert- 
^vq . Leonidov-Kovner (JIMWLK) [TTrfIB"] evolution equations were found in [T9rl22] . The question 
! of the running coupling scale in gluon production was first addressed in [23], with the corrections 
to the lowest-order gluon production cross section calculated in [21]. 

Calculations of running coupling corrections for these small- a; observables serve two related 
purposes. One is to answer the theoretical question about how running coupling effects enter 
small- x evolution and particle production. The other reason for calculating running coupling 
corrections is phenomenological: it is known that running coupling effects reduce the effective 
value of the BFKL pomeron intercept j5j[22l|25] as compared to the fixed-coupling results, 
leading to a slower growth with energy of cross sections and structure functions, which is more 
in line with the experimental data. This feature resulted in a very successful phenomenology 
developed over the past few years for the deep inelastic scattering (DIS) [26l[27] and for heavy 
ion collisions [28] based on the running coupling BK (rcBK) equation. 

The existing phenomenological applications of the running coupling evolution to DIS mainly 
concentrate on the total DIS cross section and the structure functions [SHIEZ]- A related DIS 
process is the single diffractive dissociation, in which the proton or nuclear target remains 
intact in the interaction, generating an adjacent rapidity gap, with the incoming virtual photon 
splitting into a qq pair along with several more gluons and quarks, which later fragment into 
hadrons with the net invariant mass M\. The process is illustrated in Fig. [TJ At very high 
energy / small- x the cross section of single diffractive dissociation is described by the non-linear 
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evolution equation derived in [29] in the leading-ln 1/x and large- N c approximations with fixed 
coupling. The equation was generalized beyond the large- N c limit in [3"0Tf32j. At the same time, 
until now the running coupling corrections to the equation of [29] have not been constructed. 
Similar to the running coupling corrections for other small- x observables, these corrections are 
needed first of all on the purely theoretical grounds, to determine the scale of the running 
coupling constant in the kernel of the equation of [29J. Second of all, they may be useful 
phenomenologically to analyze single diffractive dissociation in the same framework as the total 
DIS cross sections. (Note though that since the invariant mass M\ of the hadrons produced 
in the single diffractive dissociation at HERA is not very large, in practical phenomenological 
applications one often limits the virtual photon breakup to the qq pair and at most one gluon 
[331435] . thus alleviating the need for an evolution equation. We will proceed here expecting 
that such approximation will require small- x evolution corrections at the future DIS machines.) 
The goal of this work is to find the running coupling corrections to the evolution equation for 
diffractive dissociation constructed in |29|. 




\J 
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Figure 1: Single diffractive dissociation in high energy DIS. 

This letter is structured as follows. In Sec. [2] we re-derive the non-linear evolution for the 
single diffractive dissociation at fixed-coupling using a method different from the original one 
used in [29J, more similar to the one used in [32]. This new method makes inclusion of the 
running coupling corrections straightforward, which is accomplished in Sec. [3j Our final result 
is shown in (12"T1) . which, along with the kernels (Tl6|) and (ITS]) demonstrates that the running 
coupling corrections to the kernel of the non-linear evolution equation for diffraction are exactly 
the same as those for the rcBK and rcJIMWLK evolution equations. 

2 Re-deriving the fixed- coupling evolution for single diffrac- 
tive dissociation 

We begin by re-deriving the fixed-coupling evolution equation for the single diffractive dissoci- 
ation in DIS originally found in [29J. Our goal is to find the cross section crj^j for the single 
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diffractive dissociation in DIS on a nucleus A resulting in a rapidity gap and production of 
hadrons with the invariant mass M x , as shown in Fig. [TJ We denote the net rapidity interval 
between the projectile qq pair and the target nucleus by Y = In s/Q 2 , where s is the center-of- 
mass energy of the qq dipole-nucleus system and Q 2 = —q 2 is the photon virtuality. As shown 
in Fig. [JJ, if we assume that the nucleus has rapidity 0, then the projectile qq dipole would 
have rapidity Y. The rapidity gap stretches from rapidity to rapidity Yq. Assuming that 
s > M\ > Q 2 one has Y ss In s/M|. 

We are working in the kinematics when both Yq and Y —Yq are very large, such that a s Yq ~ 1 
and a s (Y — Yq) ~ 1 and the leading-logarithmic small- x evolution corrections resumming powers 
of a s Yq and a s (Y — Yq) are important. That is, we need to resum gluon emissions and virtual 
corrections for the gluons inside the rapidity gap and in the dissociation region of the qq pair 
leading to creation of hadrons with the invariant mass M\. 

Following [29] we denote by N® Xl (V, Yq) the cross section per unit impact parameter of 
the single diffractive dissociation in the dipole-nucleus scattering with the rapidity gap greater 
than or equal to [0, Yq] and the net rapidity interval of Y. As usual xo and xi are the transverse 
positions of the quark and the anti-quark in the qq dipole, as shown in Fig. [TJ (Boldface notation 
denotes two-component vectors in the transverse plane.) With the help of N D the differential 
cross section for the single diffractive dissociation in DIS can be written (for a fixed s) as 



M l^§f = - j d 2 x d 2 Xl J dz |^^(xoi,^)| 



*x J 9Y 



with {^'^^(xqi, z)\ 2 the standard order-a^A/ light-cone wave function squared for a virtual 
photon fluctuating into a qq pair with xoi = |x — xi| (see e.g. [9]). 

We want to write down an evolution equation for iV^ x (Y,Yq). First we note that when 
Yq — Y the rapidity gap spans the whole rapidity interval and the interaction is elastic. In such 
case [29] 



< )X1 (^ = Yq,Y ) = [iV X0iXl (r )] 2 (2) 
where A^" XOiXl (F) is the dipole-nucleus forward scattering amplitude defined by 

N m (Y) = j^- (tr[l-V X0 Vl]) Y (3) 

with Vx the Wilson line describing the interaction of a quark at x moving along the 
with the target 

V x = Pexp{ig I dx + t a A a ~(x + , x~ = 0, x) J> . (4) 



Here N 1 is the operator of the gluon field of target, t a are fundamental generators of SU(iV c ), 
and (. . )y in Eq. (|3J) denotes averaging in the target wave function evolved up to rapidity Y (see 
for reviews of this notation). The light cone coordinates are defined clS X — (t±z)/V2. 
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In the large- N c limit and in the leading- In s approximation A r XQjXl (Y) obeys the BK evolution 
equation [5THU] 

a N f x 2 

d Y N X(hXl (Y) = / d 2 x 2 [iV XOiX2 (F) + N X2 , K1 (Y) - iV XOiXl (F) 

Z 71 J Xq2 %2i 

-N X0lX2 (Y)N X2iX1 (Y)} (5) 

where Xij = |xy| with Xy = x, — Xj and <9y = d/dY. The solution of Eq. ((Sj) would allow us to 
find the elastic dipole-nucleus scattering cross section (per unit impact parameter) via Eq. (j2J) • 
The resulting N% (y = Y ,Y ) would serve as the initial condition for the ^-evolution of 
N x Xl (Y, Yq) we are about to construct. 

Note that N D is a scattering cross section (per unit impact parameter), and, therefore, 
when written as amplitude squared, contains interactions to the left and to the right of the 
final state cut. (For example, in Eq. (T5]) the interactions are contained in the two factors of N.) 
It is easier to derive its evolution by completing it to the S'-matrix like object by defining 

S^(Y,Y ) = l-2N X(hXl (Y) + N^ Xi (Y,Y ). (6) 

The new quantity S D contains the no-interaction term (= 1), and two terms with the interac- 
tions only to the left or to the right of the cut (giving — N XQjX1 (Y) each), in addition to N D . 
The addition of these extra terms which have no new particles produced in the final state to 
N D means that all terms in S D describe the final state with the rapidity gap greater than or 
equal to Yq, allowing for no-interaction contributions to the right and to the left of the cut. 
The quantity S D is illustrated in Fig. [2, where the solid vertical straight line denotes the final 
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Figure 2: Diagrammatic representation of S D defined in Eq. fl6]). 

state cut at the light-cone time x + = oo, while the dashed vertical straight lines denote the 
Glauber-Mueller interaction with the target [38j which happens at the time x + = both in 
the amplitude and in the complex conjugate amplitude. The times are labeled in the leftmost 
diagram in Fig. [2j (The time scale of the interaction with the target is much shorter than the 
time scales of the small-x evolution: see [36,33,139] for reviews of similar calculations.) Shaded 
ovals and the circle in Fig. [2] represent all gluon emissions and interactions with the target 
nucleus contained in S D , N, and iV D as labeled. 

To construct one step of small-x evolution for S D we have to see how this quantity changes 
under an emission of a gluon with rapidity Y > Y . Note that for Y = Y Eq. (j2j) leads to 

SZ^P = Y o, Y ) = [1 - iV Xo , Xl (F )] 2 , (7) 

which is the initial condition for the evolution of S D . In one step of small-x evolution a soft 
gluon can be emitted and absorbed by the quark and the anti-quark in the dipole 01 at any time 
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x + G (— oo, +00) both in the amplitude and in the complex conjugate amplitude. The rapidity 
gap constraint on the final state of S D implies that no gluon should be present at x + = 00 for 
gluon rapidities less than Y . However, this condition is already satisfied by our initial condition 
(JTJ). There are no final-state restrictions on the gluons with Y > Yq: they may or may not go 
through the cut. This makes the gluon emissions and absorptions at times x + G (0, +00) in the 
amplitude and in the complex conjugate amplitude cancel, due to the final-state cancellations 
originally derived in [ID] (see e.g. j29HH] for applications of such cancellations). The final state 
cancellations of [10] are illustrated in Fig. [3] which shows how the diagrams with the gluon 
either emitted or absorbed (or both) at x + G (0, +00) cancel. (Cancellations also take place 
for all other couplings of the gluon to the dipole not shown explicitly in Fig. [3j) 




Figure 3: Diagrammatic examples of the final state cancellations derived in [ID]. Cancellation 
for the sum of the diagrams which are complex conjugates of those shown in the bottom panel 
also takes place, but is not shown explicitly. 

This leaves us with gluon emissions/absorptions at x + G (— 00, 0) both in the amplitude 
and in the complex conjugate amplitude. Identifying the time interval x + G (—00, 0) in the 
complex conjugate amplitude with the time interval x + G (0, +00) in the forward amplitude 
we conclude that the evolution for S D with Y > Yq is the same as for the S'-matrix of the 
dipole-nucleus scattering, which is related to the forward amplitude N by 

S X0 , Xl (F) = l-iV X0 , Xl (r). (8) 

In the large- N c limit it is described by the BK equation (jSj) re-written for the S'-matrix: 

dvS^ (Y, Y ) = ^ / d 2 x 2 [ S ^ 2 (Y, Y ) < Xl (Y Y ) - < Xl (F, Y )] . (9) 

Z 71 J X 2 q X 2 i 

The evolution equation (J9j is illustrated in Fig. HI where, in the large- N c approximation notation 
gluons are denoted by the double lines and disconnected gluon (double) line implies all possible 
connections to the quark and the anti-quark in the dipole (see [121 - PH] ). 

One could stop at Eq. ([9]), since its solution with the initial condition flTJ) used in Eq. <^ 
would give one N D . (Even more, Eq. (J6j implies that dy S D = dY N D , which is all that is 
needed to find the differential diffractive cross section in Eq. (JTJ).) Alternatively one can use 
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Figure 4: The evolution equation for S D . 



Eq. ([6]) along with Eq. (j5j) in Eq. to obtain the evolution equation for N D [29] 

a s N c 



2tt 5 

tD 



d xi 



<-y> 2 ™ 2 

x 20 x 21 



<, X2 (v, n) + <, X1 (y, y ) - < >X1 (y y 



+<, X2 (K, y ) < iXi (y, r ) - 2 iv X0|X2 (y ) < >X1 (y, y ) - 2 < >X2 (y y ) iv X2iX1 (y ) 

+ 2iv X0iX2 (y)iv X2 , Xl (y) 



(10) 



3 Running coupling corrections 

Similar to the previous works on the running coupling corrections to small-x observables [19}- 
[23] we use the Brodsky-Lepage-Mackenzie (BLM) prescription [45] : one has to dress all the 
gluon propagators and multi-gluon vertices with quark bubble corrections, and then replace the 
number of quark flavors Nf by 



where 



02 



llNr - 2N 



f 



12 7T 



(11) 



is the one-loop QCD beta-function. In the end one is able to absorb all such corrections into 
the powers of the physical running couplings 



a s (Q 2 



a,. 



1 + a„ 2 In ^ ' 



(13) 



This procedure fixes the scales Q in the running couplings. 

When the gluon lines in our above derivation of Eq. ( TTOT) are dressed by the quark loop 
corrections, the final state cancellations of [30] shown in Fig. [3] still take place. All the g — >■ qq 
splitting and the qq — > g mergers happening at the times x + G (0, +oo) in the amplitude 
and in the complex conjugate amplitude cancel (simultaneously with the cancellation of gluon 
emissions and absorptions by the original dipole at x + G (0, +oo) shown in Fig. [3]). These 
cancellations have been discussed in detail in [23]. An example of a cancellation of the final 
state interactions with a quark loop correction added to the top line of Fig. [3] is shown in Fig. |5j 
as we demonstrate in Appendix A, the sum of all the diagrams shown in Fig. [5] is zero, 



A + B + C + D + E = 0. 



(14) 
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Figure 5: An example of the final-state interactions with a quark loop correction: the sum of 
these diagrams is zero, as shown in Appendix A. 



The argument can be repeated at higher orders, showing that higher numbers of quark loops 
would also cancel in the final state, and for other diagrams with the final-state interactions. 

The remaining g — )■ qq splitting and qq — » g mergers are limited to the x + G (— oo, 0) times 
both in the amplitude and in the complex conjugate amplitude, thus dressing the propagators 
of the early-time gluons only. Therefore, the running coupling evolution for S D is the same as 
the running coupling evolution for the BK equation (JHJ) written for the dipole-nucleus S'-matrix 
of Eq. (jSJ). Using the results of the derivation of the rcBK evolution in [191120] we write the 
running coupling evolution equation for S D as 

dySg^Yo) = J d 2 x 2 K{^,^ 2 ) K >X2 (y,Y ) < Xl (Y,y ) - S° >X1 (Y,Y )] . (15) 



The kernel iT(x , xi, x 2 ) of the evolution equation (|T5|) depends on the prescription used to fix 
the scales of the running couplings (see [22] for details). In the Balitsky prescription |20j the 
kernel is 



Af c a '(x , Xl ,x 2 ) 



N c a s (x 2 10 ) 
2tt 2 



x 
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x 20 x 21 



+ 



a s( x 2 0J 



•^20 \ a s( x 2l) 



X 21 V^s^ld) 



(16) 



with the following shorthand notation 



a s (x 2 ± ) = a s 



4 e -|-2 7B 



;i7) 



in the MS renormalization scheme. In the Kovchegov-Weigert prescription [19] the kernel is 



^(xo.xx.xa) 



It: 2 



a s {x 2 2Q )- 



s( X 2l) X 20 ' x 21 



X 



20 



a s (R 2 



^y»2 ry, 2 

x 20 x 21 



+ a s (x 



X 



21 



;is) 



with 
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X20 %21 



X21 
X20 



2,2 22 

2" _ 2 2 - 

c 20 ^21 



(19) 



The terms in the kernel neglected by each prescription can be found in [20|f22] . 
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The initial condition for the running coupling evolution for S D is given by Eq. (J7J), where 
N is now found from the rcBK equation [T91I2T)] 

d Y N X(hXl (Y) = J rf 2 x 2 ^(x , Xl ,x 2 ) [N X(hX2 (Y) + N m (Y) - N X{hXl (Y) 

-N X0 , X2 (Y)N X2iX1 (Y)]. (20) 
Finally, using Eq. (|6]) in Eq. (JT5l) one obtains the running coupling evolution for N D 

d Y N^ xi (Y,Y ) = j d 2 x 2 if(xo,x 1)X2 ) < >X2 (y,y )+iv^ xi (y ) Fo)-<, xl (^n) 

+< iX2 (Y, Y ) iV x D 2>xi (Y, Y ) - 2 iV X0 , X2 (Y) iV^ (Y, Y ) - 2 < jX2 (Y, Y ) iV X2iXl ( Y) 

+ 2iV X0 , X2 (Y)iV X2 , Xl (Y) (21) 

with the kernel given by Eq. (|T6l) or by Eq. ( 1T81) depending on the prescription choice and with 
the initial condition given by Eq. ([2]) with N provided by Eq. ( 1201) . 

Eqs. f|T5|) and fl2Tj) are the main results of this letter, providing running coupling corrections 
to the non- linear small- x evolution for the cross section of diffractive dissociation in DIS. The 
kernels in Eqs. (fT6j) and f|T8|) have of course been found before in the calculation of the rcBK 
and rcJIMWLK evolutions [19l[20]: we showed that these same kernels govern the evolution 
for diffractive dissociation. Our derivation was done in the large- N c limit: since the evolution 
equation ( |T5l) is simply the rcBK equation, its generalization to all-iV c is the same as for the 
rcBK equation, and is accomplished by performing the (. . .)y,y averaging to the whole right- 
hand side of Eq. (|T5l) instead of applying it to each individual S D separately (see [T5l[T6] ). 
Note that the averaging should now include the final-state rapidity gap greater than or equal 
to[0,Y ]. 
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A Final state cancellations 

We want to demonstrate the cancellation of diagrams in Fig. [51 We will be working in the 
light-cone perturbation theory (LCPT) [161. H7]. In LCPT the diagrams in Fig. [5] differ only by 
their energy denominators and overall signs. Suppressing the rest of the diagrams we write for 
the contribution of the diagram A 



(E, + E 2 y {Ei + E 2 - E k f 
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where Ej, = k 2 / (2 k + ) is the light cone energy of a line with momentum k, and E\ = k 2 /(2 k± ), 
£^2 — k 2 ./^ /c^~) (momentum labels are explained in Fig. 

The remaining diagrams have to be treated more carefully, due to the presence of interme- 
diate states with zero light-cone energy denominators [231H0]. Noting that the diagram B has 
an overall minus sign compared to A due to moving a quark-gluon vertex across the cut we 
write 

OO tg t2 OO 

B = — J dt 3 e i ( E k-Ei-E 2 +ie)t 3 J dt2e i(E 1+ E 2 -E k +ie)t 2 j ^ j {E k +i e) H J ^ g -i t* 



+ 0(6), (A2) 





2 E\ + 2 E2 — 3 E k 



2eE 2 k (E 1 + E 2 - E k ) 2E\ (E 1 + E 2 - E k 

where all the light cone times associated with the vertices, ti = xf, are labeled in Fig. [5j3 and 
e is an infinitesimal regulator. (See [221110] for similar calculations.) One can show that the 
singular 0(l/e) term in Eq. (IA2I) is canceled by the diagrams with the gluon lines replaced by 
the instantaneous interactions of LCPT: we thus concentrate on the 0(e°) term. Also, due to 
mirror symmetry, C = B*, such that 0(l/e) terms would cancel in the B + C sum [40J. 
Similarly for the diagrams D and E we get (see Fig. 03 for notations) 

t4 <3 t2 



D = j dUe i( ~ Ek+ie)U j dtze i{E k -E^E 2+ ie)t s J dhe i { E 1+ E2-E k+l e)t 2 J^^^ 



+ 0(e). (A3) 



000 
2E 1 + 2E 2 + E k 



AeE 2 k (E 1 + E 2 ) 2El(E 1 + E 2 ) 

Again, due to mirror symmetry, E = D*. With the help of Eqs. ( 1A1I) . ( 1A2I) . and ( 1A3I) and 
dropping the 0(e) terms we see that 

^ + 5 + + ^ + ^ = 0, (A4) 

as desired. 
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